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^ . Abstract 



The azimuthal dependence of 3D HBT radii relative to the event plane gives us information 

about the source shape at freeze-out. It also provides information on the system's evolution by 

comparing it to the initial source shape. In recent studies, higher harmonic event planes and flow 

^D I have been measured at RHIC and the LHC, which result primarily from spatial fluctuations of 

the initial density across the collision area. If the shape caused by initial fluctuations still exists at 

freeze-out, the HBT measurement relative to higher order event plane may show these features. 

^ ' We present recent results of azimuthal HBT measurements relative to 2'"'- and 3''''-order event 

1* I planes in Au-hAu 200 GeV collisions with the PHENIX experiment. Recent HBT measurements 

r— H . at lower energies will be also shown and compared with the 200 GeV result. 

1. Introduction 

\^ , HBT measurements provide information on the space-time evolution of the particle emitting 

t^^ ' source in relativistic heavy ion collisions. The azimuthal dependence of 3D HBT radii with 

OO , respect to an event plane gives us information on the source shape at freeze-out. It also provides 

^N ■ information on the system's evolution by comparing it to the initial source shape. The higher 

t:J" , harmonic flow (v3, v'4, etc) of particles has recently been measured at RHIC and the LHC. It is 

CD • primarily due to the spatial fluctuation of the initial density of the collision area. A hydrodynamic 

model calculation |1] reports that the shape of the initial fluctuations resulting in a triangular 
component of the initial shape may be preserved until freeze-out. HBT measurements relative to 
a higher-order event plane may reveal this feature. 



H ' 2. Azimuthal HBT witli respect to tlie 2"''-order event plane 



Azimuthal HBT radii with respect to the 2""'-order event plane have been measured for 
charged pions and kaons in ^/s^ = 200 GeV AuH- Au collisions at PHENIX ||2[] . It was found that 
the final eccentricity of kaons, which is defined as S2jinai - 2^^2/^iO 1211' i^ larger than that of 
pions and almost the same as the initial eccentricity. However, since HBT radii show a transverse 
mass (mj-) dependence and the average nij of pions and kaons are different, the nij dependence 
needs to be considered in the comparison of the final eccentricities. Figure [T] shows the relative 
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amplitude of the azimuthal HBT radii for charged pions and kaons as a function of {mj) for two 
centraHty bins. The left top panel corresponds to the final eccentricity; there is still a difference 
between pions and kaons in non-central collisions even at the same {mr). This difference may 
be due to different cross-sections of pions and kaons. The relative amplitudes of R„ and Ros at 
low nij in the most central collisions have finite values although the final eccentricity is almost 
close to zero. This result may indicate a temporal variation of the emission duration of particles 
because Rg and Ros contain temporal information in addition to geometrical information. 
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Figure 1 : Relative amplitude of azimuthal HBT radii for charged pions and kaons with respect to 2 -order event plane 
in Au-l-Au 200 GeV colUsions. 



3. Azimuthal HBT with respect to 3'^''-order event plane 



Fluctuations in the initial geometry of the heavy ion collisions, which is considered to be the 
origin of higher harmonic flow, may be preserved until freeze -out. Triangular flow, V3, is known 
to have a weak centrality dependence, while the initial triangularity calculated within a Glauber 
model has a pronounced centrality dependence [4]. Triangularity at freeze-out is determined by 
the initial triangularity, V3, the expansion time, and so on, and will provide detailed information 
on the space-time evolution. 

Figure|2]shows R^ and R,, for charged pions as a function of the azimuthal angle with respect 
to 2'"'- and 3''''-order event planes (^2 and T3) in 0-10% in Au-t-Au 200 GeV collisions, where 
the averages of the radii with respect to T2 and T3 are set to 10 and 5 fm^ respectively. Filled 
symbols show measured data points, and open symbols are the same data points reflected around 
<p - ^'„ = 0. The solid lines depict the fit functions 7?^_o+2/?/j,«cos(n(0 - ^„)) {jj. - s, o). The 
values of R^ show a very weak oscillation with respect to both ^2 and *I'3, while Rg clearly 
exhibits stronger oscillation. The oscillation strength of R„ relative to *P3 is comparable to that 
relative to "^2- The oscillation of Rg may carry information about the duration of the emission. 
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Figure 2: The azimuthal dependence of R, and /?„ for charged pions with respect to 2'"'- and 3"'-order event plane in 
Au-hAu 200 GeV collisions, where the averages of radii with respect to 'Ft and "F3 are set to 10 and 5 fm^ respectively. 



The ^2 dependence may indicate that the emission duration is different in-plane versus out-of- 
plane since R,, shows a weak oscillation and the source shape is thought to be close to a circle. 
The *P3 dependence may be also due to the difference of emission duration between different 
azimuthal directions. However, the flatness of Rs doesn't necessarily mean that the source shape 
is not triangular but circular ylj. It may be difficult to imagine that the emission duration has 
such strong variations in azimuth relative to different event planes. The oscillation of Rg may 
reflect not only the emission duration but also the depth of the elliptical or triangular source. 

As mentioned before, the final eccentricity is defined as S2ji„ai - 2^j 2/^10- Here, we define 
B-ijinai as IR^^T^jR^^^. Figure |3] shows s„jinai as a function of e„_,„,„a;, where s^jninai is calculated 
using a Glauber model. Note that S3jmai doesn't represent the final triangularity because there 
is no higher harmonic anisotropy in the Gaussian approximation for a static source. However, 
it will mean any triangular component of homogeneity region in a expanding source. The ob- 
served E-ijinai doesn't seem to exhibit any centrality dependence and is zero within systematic 
uncertainties. 
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Figure 3: Final £„ as a function of initial e„ in Au4-Au 200 GeV collisions. Initial £„ is calculated in a Glauber model. 
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Figure 4: The product of pion 3D HBT radii as a function of charged particle multiplicity density. 

4. Low energies at PHENIX 

The RHIC beam energy scan program has been conducted in order to explore the critical 
point between the QGP and hadron gas phases in the QCD phase diagram, The centrality and 
niT dependence of HBT radii were measured for 39, 62 and 200 GeV collision energies. We 
observe no significant change beyond systematic errors across the three energies. Figure |4]shows 
the product of 3D HBT radii for charged pions as a function of charged multiplicity density. 
PHENIX results are compared with results at different energies and collision species at AGSl|5|], 
SPS||6|,i3l, RHIcfs','?] and LHCfl^. The product of 3D HBT radii, which represents the volume 
of the homogeneity region, from PHENIX is consistent with the global trend. 

5. Summary 

The latest results of azimuthal HBT measurements with respect to the 2"''- and 3"'-order event 
planes are presented. The difference of final eccentricity for charged pions and kaons is seen 
even in the same mj- region, which may indicate a shorter freeze-out time of kaons due to a lower 
cross section. The azimuthal dependence of pion HBT radii relative to *I'3 has been measured 
and the oscillation of R,, is clearly seen. This may be indicative of the temporal variation of the 
emission duration or the depth of the source with triangular shape. These results will provide new 
constraints on theoretical models and new insight into the space-time evolution of the system. 
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